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ABSTRACT 

We examine whether or not it is possible to derive the field strength distribution of quiet Sun 
internetwork regions from very high spatial resolution polar imetric observations in the visible. In 
particular, we consider the case of the spectropolarimeter attached to the Solar Optical Telescope 
aboard Hinode. Radiative magneto-convection simulations are used to synthesize the four Stokes 
profiles of the Fe I 630.2 nm lines. Once the profiles are degraded to a spatial resolution of ff!32 
and added noise, we infer the atmospheric parameters by means of Milne-Eddington inversions. The 
comparison of the derived values with the real ones indicates that the visible lines yield correct 
internetwork field strengths and magnetic fluxes, with uncertainties smaller than ~150 G, when a 
stray light contamination factor is included in the inversion. Contrary to the results of ground-based 
observations at f , weak fields are retrieved wherever the field is weak in the simulation. 
Subject headings: Sun: magnetic fields - Sun: photosphere - Instrumentation: high angular resolution 



4^ . 1. INTRODUCTION 

I The characterization of quiet sun internetwork (IN) 

O fields is an important issue in solar physics. Polar imetric 
|H I measurements of the visible Fe I lines at 630.2 nm and 
^ the near-infrared Fe I lines at 1565 nm have been used 
, to advance our knowledge of IN magnetism, but no con- 

sensus has been reached yet. While the analysis of the 
visible lines suggest a predominance of kC field strengths 
^ ■ and small filling factors (Sanchez Almeida & Lites 2000; 
^sD ' Dominguez Cerdeha et al. 2003; Socas-Navarro & Lites 
On 2004), the near-infrared lines indicate that most fields 
' have hC strengths with larger filling factors (Lin 1995; 
O ; Lin & Rimele 1999; Khomenko et al. 2003; Martinez 
l/^ . Gonzalez et al. 2006; Dominguez Cerdeha et al. 2006). 
■ Attempts to reconcile these contradictory results have 
argued that visible and IR lines sample different mag- 
netic structures in the resolution element (Sanchez 
Almeida & Lites 2000; Socas-Navarro & Sanchez 
. ^ Almeida 2003) or that noise affects the visible lines more 
dramatically than the IR lines (Bellot Rubio & Colla- 
dos 2003). On the other hand, Martinez Gonzalez et al. 
. 5^ (2006) have convincingly demonstrated that it is not pos- 
sible to obtain reliable IN field strengths from the Fe I 
630.2 nm lines due to crosstalk with thermodynamical 
parameters at 1" resolution. 

Different techniques have been proposed to study IN 
fields, but the diagnostic potential of high spatial reso- 
lution observations in the absence of atmospheric seeing 
has only been explored by Khomenko et al. (2004). The 
spectropolarimeter (SP; Lites et al. 2001) of the Solar 
Optical Telescope (SOT) aboard Hinode (Ichimoto et al. 
2005) is already providing nearly diffraction-limited ob- 
servations of the solar photosphere, and upcoming in- 
struments will do so in the future (e.g., IMaX aboard 
SUNRISE or VIM aboard Solar Orbiter). Thus, there is 
a clear need to assess whether reliable IN field strengths 
can be derived from space-borne observations. 

Here we address this question by simulating and an- 
alyzing Hinode/SP measurements. Radiative magneto- 
convection simulations are used to synthesize the Stokes 
profiles of the Fe I 630.2 nm lines. The profiles are 



degraded to the nearly diffraction-limited resolution of 
0'.'32 achieved by Hinode/SP. After adding noise to the 
Stokes spectra, the atmospheric parameters are inferred 
by means of inversion techniques. By comparing the de- 
rived values with the real ones we determine the errors in 
field strength and magnetic flux to be expected from the 
analysis of Hinode data. Our main result is that Milne- 
Eddington (ME) inversions of the visible 630.2 nm lines 
yield correct IN field strengths with uncertainties smaller 
than 150 G for the whole range of strengths from ~0.1 
to 1 kG. If internetwork fields are hG fields, then simple 
ME inversions of Hinode/SP measurements will result in 
sub-kG field strength distributions, contrary to what is 
obtained from current ground-based observations. 



2. 



MHD SIMULATIONS, SPECTRAL SYNTHESIS, AND 
IMAGE DEGRADATION 



To describe the sun in the more realistic way possi- 
ble we use radiative MHD simulations by Vogler et al. 
(2005). Specifically, we take three snapshots from dif- 
ferent simulation runs representing very quiet, unipolar 
internetwork and network regions with average unsigned 
fluxes of 10, 50, and 200 Mx cm~^, respectively. The 
horizontal and vertical extents of the computational box 
are 6 and 1.4 Mm. The synthesis of the Stokes spectra of 
the two Fe I lines is carried out using the SIR code (Ruiz 
Cobo & del Toro Iniesta 1992). The spectral region is 
sampled at 113 wavelength positions in steps of 2.15 pm, 
following the Hinode/SP normal map mode (for details, 
see Shimizu 2004). The atomic parameters have been 
taken from the VALD database (Piskunov et al. 1995). 

The aperture of Hinode/SOT is 0.5 m, which operating 
at 630 nm provides a spatial resolution of ~0'.'26 (equiv- 
alent to ~ 190 km on the solar surface). The sampling 
interval in the MHD simulations is 0'.'0287, implying a 
spatial resolution of 0'.'057 (41.6 km). Thus, in order 
to simulate SOT observations, the synthetic Stokes pro- 
files derived from the MHD model have been spatially 
degraded by telescope diffraction. In addition, we have 
considered the extra loss of contrast caused by the inte- 
gration of the signal in the detector. Finally, the images 
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Fig. 1. — Left and middle: Continuum intensity maps for the simulation snapshot with average unsigned flux of 10 Mx cm~^ and for 
the data spatially degraded considering telescope diffraction and pixel size. Color scales are the same in the two maps. The contrast varies 
from 13.7% in the original image to 8.5% in the degraded one. Right: MTF of the detector (dotted line), diffraction limited MTF (dashed 
line), and combination of both effects (solid line). 

gion (Westendorp Plaza et al. 2001; Bellot Rubio 2006). 

We apply the ME inversion to the Fe I 630.15 nm and 
Fe I 630.25 nm lines simultaneously. A total of 9 free pa- 
rameters are determined (So, Si, t^o, AA^, a, 5, 7, x, and 
vlos; for the meaning of the symbols see, e.g., Orozco 
Suarez & del Toro Iniesta 2007). No additional broad- 
ening of the profiles by macroturbulence or microturbu- 
lence is allowed. Three different inversions are performed 
to derive the atmospheric parameters. All of them use 
a simple one-component model, i.e., a laterally homoge- 
neous magnetic atmosphere occupying the whole resolu- 
tion element. We first invert the profiles in the absence of 
noise, and then with noise added at the level of 10~^/c- 
In the last inversion, the noisy profiles are fitted consid- 
ering non-zero stray light contaminations factors. The 
stray light profile is evaluated individually for each pixel 
by averaging the Stokes / profiles within a box l^'-wide 
centered on the pixel. For all inversions we use the same 
initial guess model, allowing a maximum of 300 itera- 
tions. The initial field strength is 100 G. 

4. RESULTS 

Figure [3] shows the vector magnetic field (strength, 
inclination, and azimuth) retrieved from the inversions 
of the Stokes profiles. The first column displays a cut 
of the simulation snapshot with average flux density of 
10 Mx cm~^ at optical depth logr = —2. The second 
and third columns contain the results of the ME inver- 
sions of the spatially degraded profiles in the absence 
of noise and the specific case of a SNR of 1000, respec- 
tively. Finally, the fourth column shows the atmospheric 
parameters derived from the noisy profiles accounting for 
stray light contamination. White regions represent pixels 
which have not been inverted because of their small po- 
larization signals (we only consider pixels whose Stokes 
Q^U otV amplitudes exceed three times the noise level). 

Over the granules, the magnetic field strength is very 
weak and the polarization signals are buried in the noise. 
These pixels represent ~55% of the total area (white re- 
gions in Fig.[3|). The stronger fields concentrate in inter- 
granular regions. In those regions, the magnetic struc- 
tures inferred from the inversion have bigger sizes than 
the real ones, i.e., they appear "blurred". This is caused 
by the degradation of the images due to telescope diffrac- 
tion and CCD pixel size. The field inclination and az- 
imuth structures resulting from the inversion are blurred 
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Fig. 2. — Comparison between the average Stokes / profile from 
the spatially degraded data (solid) and the FTS spectral atlas 
(dashed). Both continua are normalized to unity. 

are rebinned to the SP CCD pixel size of 0Vl6x 0'.'16. 
Figure [T] shows continuum intensity maps for the orig- 
inal and the spatially degraded data. The degradation 
process reduces the rms contrast from ~14% to ~9% in 
the continuum. The pixelation of the CCD is notice- 
able in the degraded image. The right panel of Fig. [T] 
shows Modulation Transfer Functions (MTFs) describing 
the filtering of spectral components induced by telescope 
diffraction and pixelation effects in the CCD. Note the 
modification of the effective MTF caused by the central 
obscuration of the entrance pupil. 

Figure [2] compares the average Stokes / profiles from 
the spatially degraded data and the NSO Fourier Trans- 
form Spectrometer Atlas of the quiet Sun. Both spectra 
are very similar, with only small differences in the line 
core and wings of Fe I 630.1 nm. Note that the lack of a 
temporal average exclude, for instance, the effect of the 
5-min oscillation in the simulated profile. To reproduce 
an Hinode/SP observation we have also convolved the 
profiles with a Gaussian of 30 mA FWHM to account 
for the spectral resolving power of the spectrograph and 
have added noise at the level of 10~^ of the continuum 
intensity Ic (the polarimetric sensitivity of standard Hin- 
ode/SP measurements). 

3. INVERSION 

To derive the magnetic field strength from the simu- 
lated profiles we use a least-square inversion technique 
based on ME atmospheres. ME inversions represent the 
best option to interpret the measurements if one is not 
interested in vertical gradients of the physical quantities. 
They are simple and often provide reasonable averages of 
the atmospheric parameters over the line formation re- 
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Fig. 3. — Left: Cuts at optical depth logr = —2 of the model atmospheres provided by the MHD simulation with averaged unsigned 
flux of 10 Mx cm~^. Second column: Maps of the physical quantities retrieved from the ME inversion of the proflles with no noise. Third 
column: Maps retrieved from the ME inversion of the proflles with SNR of 1000 and no stray light contamination. Fourth column: Same 
as before, but accounting for stray light. From top to bottom: magnetic fleld strength, inclination, and azimuth. 



as well. The azimuth values are rather uncertain be- 
cause of the tiny linear polarization signals produced by 
the weak fields of the simulations. 

Figure [4] is a close up of small features observed in 
intergranular regions. Note that each Hinode/SP pixel 
of 0V16 X 0V16 corresponds to 36 pixels in the simula- 
tion, hence they usually contain a broad distribution of 
magnetic field strengths. When we consider that the po- 
larization signal is produced by a single magnetic compo- 
nent within the resolution element and no stray light is 
allowed for, the inferred field strengths are smaller than 
those in the model, so the field is underestimated (middle 
panels of Fig. 2]). If one accounts for stray light contami- 
nation the inferred fields become stronger (right panels), 
but also noisier due to the increased number of free pa- 
rameters. 

To analyze the results in a more quantitative way we 
calculate the mean and rms values of the errors. We de- 
fine the error as the difference between the inferred and 
the real parameters at optical depth logr = — 2. Since 
one pixel of the degraded data corresponds to 36 pixels 
in the simulations, we compare each inverted pixel with 
the mean of the corresponding 36 pixels in the original 
map. Figure [5] shows the mean and rms errors of the 
field strength resulting from the inversion without ac- 
counting for stray light (top left panel). It is clear that 
fields above ~100 G are underestimated, with rms errors 
smaller than ~150 G in the whole range of strengths. 
The results are similar for the magnetic fiux density (top 
right panel). The inversion considering stray light con- 
tamination yields much better inferences, as can be seen 
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Fig. 4. — Left: Field strengths at logr = -2 in the MHD simula- 
tions with 10 Mx cm~^ (top) and 50 Mx cm~^ (bottom). Middle: 
Field strengths derived from the ME inversion of the spatially de- 
graded Stokes proflles with SNR 1000 and no stray light contam- 
ination. Right: Field strengths from the ME inversion accounting 
for stray light contamination. 

in the bottom panels of Fig. [5l The field strength and 
fiux are slightly overestimated for weak fields, but the 
rms errors do not exceed 150 G in any case. 

5. DISCUSSION 

Previous analyses of visible (630.2 nm) and near- 
infrared (1565 nm) neutral iron lines do not agree on 
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Fig. 5. — Top: Mean (solid) and rms (dashed) errors of the field 
strength (left) and flux (right) derived from the inversion of the 
proflles with SNR=1000 assuming a single magnetic atmosphere 
and no stray light contamination. Bottom: Same as before but 
accounting for stray light contamination. 

the distribution of field strengtiis in IN regions. In par- 
ticular, the visible lines systematically deliver kG field 
strengths and small filling factors, while the near-infrared 
lines suggest a predominance of hG fields. 

Here we have shown that ME inversions of the 
Fe I 630 nm lines at spatial resolutions of 0^32 (the case 
of Hinode/SP) underestimate the magnetic field strength 
by some hundred G if no stray light contamination is in- 
cluded. When stray light is accounted for, ME inversions 
are able to recover any magnetic field above 100 G with 
remarkable accuracy. 

Interestingly, we always derive weak fields from the 
simulated Hinode/SP observations where the field in the 
MHD model is weak. Likewise, pixels assigned strong 
fields by the ME inversion actually correspond to strong 
fields in the MHD model. This is in sharp contrast with 
the results of Martinez Gonzalez et al. (2006), who always 
infer kG fields from the Fe I 630.2 nm lines observed in 
the IN at resolutions of 1 5 when the inversion is ini- 
tialized with strong fields. The difference is probably due 
to: (a) our significantly higher spatial resolution, which 
narrows the range of field strengths present in the pixel; 
(b) the fact that we do not employ two-component at- 
mospheres, micro- or macro-turbulent velocities, which 



reduces the degrees of freedom of the solution; and (c) 
the simple description of the thermodynamics provided 
by the ME model which, contrary to the atmosphere 
used by Martinez Gonzalez et al. (2006), does not allow 
to compensate for incorrect magnetic parameters. On 
the other hand, our results seem to contradict the con- 
clusions of Bellot Rubio & Collados (2003). However, 
the signals considered here are larger by a factor of ^10 
due to the much higher angular resolution (which implies 
larger filling factors). Under these conditions, noise does 
not significantly affect the field strengths derived from 
the visible lines (see Fig. 5 in Bellot Rubio & Collados 
2003). 

We caution that the results of this Letter may only be 
valid as long as the MHD simulations provide a realistic 
description of the Sun. The performance of ME inver- 
sions could be different if the magnetic field is structured 
on scales much smaller than ~0V3. For the moment, how- 
ever, there is no compelling evidence that tiny magnetic 
elements exist in the quiet solar photosphere. 

6. CONCLUSIONS 

Our analysis suggests that Hinode/SP observations 
will make it possible to determine the real distribution 
of field strengths in quiet Sun internetwork regions. Sim- 
ple one-component Milne-Eddington inversions without 
macro and microturbulence seem appropriate to achieve 
that goal. However, it will be essential to account for the 
degradation of the image induced by telescope diffraction 
and detector pixel size. The work presented here shows 
that the effects of the degradation can be modeled suffi- 
ciently well including a stray light contamination factor 
in the inversion. 
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